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^ . Abstract 

' The renormalization-group functions of the two-dimensional n- 

vector A(/)^ model are calculated in the five-loop approximation. Per- 
ff^ I turbative series for the /3-function and critical exponents are resummed 

■ by the Pade-Borel-Leroy techniques. An account for the five-loop term 

shifts the Wilson fixed point location only briefly, leaving it outside 
the segment formed by the results of the lattice calculations. This is 
argued to reflect the influence of the non-analytical contribution to 
• the /3- function. The evaluation of the critical exponents for n = 1, 

n = and n = — 1 in the five-loop approximation and comparison 
of the results with known exact values confirm the conclusion that 
non-analytical contributions are visible in two dimensions. For the 
2D Ising model, the estimate uj = 1.31(3) for the correction-to-scaling 
exponent is found that is close to the values resulting from the high- 
temperature expansions. 

The field-theoretical renormalization-group (RG) approach proved to be 
a powerful tool for calculating the critical exponents and other universal 
quantities of the basic three-dimensional (3D) models of phase transitions. 
Today, many- loop RG expansions for /5-functions (six-loop), critical expo- 
nents (seven- loop), higher-order couplings (four-loop), etc. of the 3D 0(n)- 
symmetric and cubic models are known resulting in the high-precision numer- 
ical estimates for experimentally accessible quantities [1-7]. The main aim of 
this report is to elaborate further the field-theoretical RG technique, namely. 
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to clear up how effective is tfiis macliinery in two dimensions wliere i) tfie RG 
series are stronger divergent and ii) singular (non-analytic) contributions to 
RG functions are expected to be larger than for 3D systems. 

The Hamiltonian of the model describing the critical behaviour of various 
two-dimensional systems reads: 



(fx 



iK^^ + (v^j^) + A(^2)2 



(1) 



where is a real n- vector field, ml is proportional to T — T^^\ T^^^ being 
the mean-field transition temperature. 

We calculate the /3-function and the critical exponents for the model (1) 
within the massive theory. The Green function, the four-point vertex and 
the 0^ insertion are normalized in a conventional way: 



= 1, (2) 

p2=0 



Gr {0,m,g4) , — 

Tr{Q, 0, 0, m, g) = mV, rj^^(0, 0, m, 5(4) = 1. 

Since the four-loop RG expansions at n = 1 are known [1] we are in 
a position to find corresponding series for arbitrary n and to calculate the 
five-loop terms. The results of our calculations are as follows: 

^ ^-g + g^- 3 f 10.33501055n + 47.67505273^ 
2 ^ ^ (n + 8)n / 

+ , ^ (^5.00027592871^ + 149.1518586n + 524.3766023^ 

-— - f 0.0888429n^ + 179.6975971^ + 2611.154871 + 7591.1087^ 

(rz + 8)H / 

+— f-0.004087i^ + 80.309671^ + 52 53.5671^ + 53218.671 + 133972 , 

n + 8H J 



\ (3) 



77 + 2 

7-^ = 1 — g + ^ {n + 2) 3.375628955 



n + 8 {n + Sy 



-,3 



^ (4.66188477271^ + 34.4184832971 + 50.18942749 ) 

(n + 8 3 V / 
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+ , ^ ( 0.3189930n=^ + 71.703302n^ + 429.42449n + 574.58772^ 



(n + 8)4 
^ 



0.09381n^ + 85.4975n^ + 1812.19n^ + 8453.70n + 10341.1^ . (4) 



2 3 

rj = ^ (n + 2) 0.9170859698 - ^ (n + 2) 0.05460897758 



+ 



^——r f-0.09268446n2 + 4.056410571^ + 29.251167n + 41.535216^ 



(n+8)n ■ y 
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(n + 8)^ 



0.07092n* + 1.05240n^ + 57.7615n^ + 325.329n + 426.896 . (5) 



In more details, the calculations are described in Ref.[8]. Instead of the 
renormalized coupling constant g4, a rescaled coupling 

9 - -^9., (6) 

is used as an argument in above RG series. This variable is more convenient 
since it does not go to zero under n ^ oo but approaches the finite value 
equal to unity. 

To evaluate the Wilson fixed point location g* and numerical values of 
the critical exponents, the resummation procedure based on the Borel-Leroy 
transformation 

oo °9 oo 

fix) = Y: = J e-H'F{xt)dt, F{y) = -^f , (7) 



i=0 



is used. The analytical extension of the Borel transforms is performed by 
exploiting relevant Pade approximants [L/M]. In particular, four subsequent 
diagonal and near-diagonal approximants [1/1], [2/1], [2/2], and [3/2] turn 
out to lead to numerical estimates for g* which rapidly converge, via damped 
oscillations, to the asymptotic values; this is clcary seen from Table 1. These 
asymptotic values, i. e. the final five-loop RG estimates for g* are presented 
in Table 2 for < n < 8 (to avoid confusions, let us note that models with 
n > 2 possessing no ordered phase are studied here only as polygons for 
testing the numerical power of the perturbative RG technique) . As Table 2 
demonstrates, the numbers obtained differ appreciably from numerical esti- 
mates for g* given by the lattice and Monte Carlo calculations [9-15]; such 
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estimates arc usually extracted from the data obtained for the linear (%) and 
non- linear (X4) susceptibilities related to each another via g^: 

2 -2 

X4 = ^TTT = -X m ^4, (8) 

Since the convergence of the numerical estimates for g* given by the re- 
summcd RG series is oscillatory, an account for higher-order (six-loop, seven- 
loop, etc.) terms in the expansion (3) will not avoid this discrepancy [8]. That 
is why we believe that it reflects the influence of the singular (non-analytical) 
contribution to the /3- function. 

The critical exponents for the Ising model {n = 1) and for those with 
n = and n = —1 are estimated by the Pade-Borel summation of the five- 
loop expansions (4), (5) for 7"^ rj. Both the five-loop RG (Table 1) and the 
lattice (Table 2) estimates for g* are used in the course of the critical expo- 
nent evaluation. To get an idea about an accuracy of the numerical results 
obtained the exponents are estimated using different Fade approximants, un- 
der various values of the shift parameter b, etc. In particular, the exponent 
7] is estimates in two principally different ways: by direct summation of the 
series (5) and via the resummation of RG expansions for exponents 

7^(2) = 1 + ^-2, T]^^^ = - - 2, (9) 

which possess a regular structure favouring the rapid convergence of the 
iteration procedure. The typical error bar thus found is about 0.05. 

The results obtained are collected in Table 3. As is seen, for small ex- 
ponent 77 and in some other cases the differences between the five-loop RG 
estimates and known exact values of the critical exponents exceed the error 
bar mentioned. Moreover, in the five-loop approximation the correction-to- 
scaling exponent uo of the 2D Ising model is found to be close to the value 
4/3 predicted by the conformal theory [16] and to the estimate 1.35 ± 0.25 
extracted from the high-temperature expansions [17] but differs markedly 
from the exact value uj — 1 [18] and contradicts to recent conjecture uj — 2 
[19]. This may be considered as an argument in favour of the conclusion that 
non-analytical contributions are visible in two dimensions. 

We thank B. N. Shalaev for numerous useful discussions of the critical 
thermodynamics of 2D systems. The work was supported by the Ministry 
of Education of Russian Federation (Grant 97-14.2-16), by the International 
Science Foundation (A. I. S., Grant 99-943), and by Saint Petersburg Ad- 
ministration (E. V. O., grant ASF 298496). 
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Table 1: The Wilson fixed point coordinate for models witli n = 1, n = 
and n = — 1 in four subsequent RG approximations and the final five-loop 
estimates for g*{n). 



n 


[1/1] 


[2/1] 


[2/2] 


[3/2] 


g*, 5-loop 


1 


2.4246 


1.7508 


1.8453 


1.8286 


1.837 ± 0.03 





2.5431 


1.7587 


1.8743 


1.8402 


1.86 ± 0.04 


-1 


2.6178 


1.7353 


1.8758 


1.8278 


1.85 ± 0.05 



Table 2: The Wilson fixed point coordinate g* and critical exponent uj for 
< n < 8 obtained in the five-loop RG approximation. The values of g* 
extracted from high-temperature (HT) [10,12] and strong couphng (SC) [11] 
expansions, found by Monte Carlo simulations (MC) [13,15], obtained by 
the constrained resummation of the e-expansion for g* (e-exp.) [12], and 
given by corresponding 1/n-expansion (1/n-exp.) [12] are also presented for 
comparison. 



n 





1 


2 3 4 


8 


9* 


RG, 5-loop 


1.86(4) 


1.837(30) 


1.80(3) 


1.75(2) 


1.70(2) 
(6=1) 


1.52(1) 
(6 = 1) 


HT 


1.679(3) 


1.754(1) 


1.81(1) 


1.724(9) 


1.655(16) 




MC 




1.71(12) 


1.76(3) 


1.73(3) 






SC 


1.673(8) 


1.746(8) 


1.81(2) 


1.73(4) 






e-exp. 


1.69(7) 


1.75(5) 


1.79(3) 


1.72(2) 


1.64(2) 


1.45(2) 


1/n-exp. 








1.758 


1.698 


1.479 




RG, 5-loop 


1.31(3) 1.31(3) 1.32(3) 1.33(2) 1.37(3) 1.50(2) 
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Table 3: Critical exponents for n = 1, n = 0, and n = —1 obtained 
Pade-Borel summation of the five-loop RG expansions for 7"^ and 
known exact values of these exponents are presented for comparison 



via the 
f]. The 



n 




9* 


7 


V 


v 


a 


1 


RG 


1.837 


1.790 


0.146 


0.966 


0.068 






1.754 (HT) 


1.739 


0.131 


0.931 


0.139 




exact 




7/4 
(1.75) 


1/4 
(0.25) 


1 








RG 


1.86 


1.449 


0.128 


0.774 


0.452 






1.679 (HT) 


1.402 


0.101 


0.738 


0.524 




exact 




43/32 
(1.34375) 


5/24 
(0.20833) 


3/4 
(0.75) 


1/2 
(0.5) 


-1 


RG 


1.85 


1.184 


0.082 


0.617 


0.765 






1.473 (SC) 


1.155 


0.049 


0.592 


0.816 




exact 




37/32 
(1.15625) 


3/20 
(0.15) 


5/8 
(0.625) 


3/4 
(0.75) 
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